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SUMMARY 


Bie  stability  of  liquid  layers  with  a  prescribed  Interface  geometry 
In  the  presence  of  surface  tension  and  unidirectional  body  forces  is 
investigated  theoretically  by  means  of  a  small  vibration  analysis. 

15ie  fluid  of  the  layers  Is  assumed  to  be  incompressible  and  Invlscld  and 
the  flow  irrotatlonal.  Particular  emphasis  is  given  to  the  effects  of 
geometry  and  unidirectional  body  forces  on  the  stability  of  such  layers. 
They  are  studied  by  means  of  a  semi- Inverse  method  which  allows  the 
exact  determination  of  particular  eigenvibrations. 

Following  the  classical  method  of  investigating  the  stability  of 
systems,  we  establish  the  equations  of  motion  for  small  perturbations 
of  the  equilibrium  configuration  and  formulate  the  boundary  conditions 
for  both  two-  and  three-dimensional  layers  with  interface  surfaces 
having  a  constant  mean  curvature.  We  linearize  the  pertinent  equations 
and  express  them  in  terms  of  the  velocity  potential.  separating  the 
time  and  space  variables,  the  vibration  problem  is  reduced  to  an  eigen¬ 
value  problem.  Exact  eigenvibrations  are  sought  by  considering  solutions 
to  Laplace's  equation  which  contain  a  number  of  free  parameters.  The 
arbitrariness  of  these  parameters  Is  reduced  by  enforcing  the  solution 
to  satisfy  the  Interface  boundary  condition.  Integrating  the  differential 
equation  of  the  streamlines  defined  by  such  solutions  and  requiring  the 
solid  supporting  surface  to  constitute  a  streamline,  we  arrive  at  rigorous 
vibration  solutions.  If  we  Introduce  sufficiently  many  free  parameters 
we  are  in  a  position  to  review  exact  vibration  solutions  corresponding 
to  a  family  of  liquid  layers  with  the  same  interface  surface  and  differ¬ 
ent  solid  supporting  surfaces  and  thus  different  thickness  distributions. 

The  stability  criteria  of  the  layers  are  presented  in  terms  of 
inequalities  involving  the  magnitude  of  either  the  unidirectional  body 

p 

force,  g,  or  the  nondlmensional  stability  parameter  n  =  T/pga  in 
which  T  denotes  the  surface  tension,  p  the  fluid  density  and  "a" 
a  characteristic  length  of  the  liquid  layer.  It  is  shown  that  the 
variation  of  the  thickness  of  a  fluid  layer  may  significantly  Influence 
its  stability,  especially  for  body  forces  in  the  near-zero  region.  The 
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results  also  reveal  that  for  some  solid  supixorting  stirfaces,  layers  of 
the  type  considered  are  only  stable  If  the  body  forces  are  larger  than 
a  certain  positive  value,  while  for  other  surfaces  stability  prevails 
as  long  as  the  body  force  exceeds  a  certain  negative  limit  value.  This 
fact  may  be  utilized  In  the  design  of  containers  of  fluid  for  near-zero 
gravity  conditions  In  such  a  manner  that  the  fluid  adheres  to  a  certain 
area  of  the  container. 

The  interpretation  of  the  results  obtained  in  this  analysis  is 
subject  to  the  restriction  that  the  dimensions  of  the  interface  surface 
be  large  compared  with  the  effective  range  of  the  contact  forces,  since 
no  contact  angle  considerations  have  been  introduced  in  explicit  form. 
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NOTATION 


Constant  coefficients 

Coefficients  of  infinite  or  finite  seriea 
Arbitrary  function  of  time 
Fimctlons  introduced  to  separate  variables 
Mean  curvature 

H  = 

m*^  order  Bessel  function  of  first  kind 

Associated  Legendre  functions  of  degree  m  and 
order 

Principal  radii  of  curvature  of  the  interface  surface 

Bounding  surface  (interface  surface  equation  =  0; 
rigid  surface  equation  =  O) 

Surface  tension 

Radius  of  equilibrium  interface  surface  or  nondimens ional 
depth  of  fluid  layer  with  flat  interface  surface  at 
X  =  0  or  r  =  0 

Half  length  of  finite  flat  uniform  fluid  layer 
Magnitude  of  body  forces 
Characteristic  length  of  layers 

Constants  used  in  separation  of  variables 

Body  force 

Integers 

Unit  vector  normal  to  rigid  surface 
Pressure  in  fluid 

Ambient  pressure  distribution  at  equlllbrixim  Interface 
surface 
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Notation  (cont'd) 


r,0 


r,z,e 


r,e,(p 


r 


0 


t 

V 

x,z 

x,z 


« 

a 

a 

a 


p 

a 


Polar  coordinates 
Cylindrical  coordinates 
Spherical  coordinates 

Nondimens  Iona  1  radius  of  axialliy  symmetric  flat  inter¬ 
face,  also  nondimens ional  radius  at  ©  =  0  of  the 
rigid  supporting  surface  Sg  of  layers  with  curved 
Interface 

Time 

Fluid  velocity 
Rectangular  coordinates 

^  X  ^  z 

Non-dimensional  length  parameters  x  =  ^  ,  ^ 

2Xq  =  non-dimensional  length  of  two-dimensional  fluid 
layers  with  flat  interface  at  z  =  0 

Velocity  potential 

Potential  function  in  terms  of  Bi>ace  variables 
Potential  of  body  forces 

Largest  angle  subtended  by  layers  with  curved  interface 
Coefficient  of  infinite  series 


oa 

g 


First  order  small  perturbation  of  interface  surface 
Free  parameter 

T  T 

Stability  parameter  pi  =  .  g  ,  or  g 

Pgh  pga 

Density  of  fluid 

Frequency  parameter  (square  of  eigenfrequency) 
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I.  INTRODUCTION 


Cie  significance  of  the  phenomenon  of  interfacial  instability  of 

liquid  layers  in  certain  ablative  heat  protection  problems,  as  well  as 

in  the  behavior  of  fluids  in  containers  under  near-zero  gravity  condi- 

( 1-3) 

tlons,  has  been  repeatedly  pointed  out.  '  The  theoretical  analysis 
of  this  phenomenon  has  so  far  been  largely  restricted  to  layers  of  ex¬ 
tremely  simple  geometry,  in  particular  to  layers  with  a  uniform 
thickness. In  this  investigation  we  consider  layers  with  a  non- 
uniform  thickness  and  interface  surfaces  with  a  constant  mean  curvature. 

With  the  aid  of  the  results  of  a  recent  study  of  the  effects  of  curvature 

(2) 

and  unidirectional  body  forces  on  the  stability  of  liquid  layers  and 

the  application  of  a  semi-inverse  method  of  flndirig  exact  solutions  to 

(12) 

the  free  oscillation  problem  of  a  fluid  in  a  container  we  derive 
stability  criteria  without  resorting  to  approximate  methods. 

As  in  reference  (2),  the  density  of  fluid  separated  by  the  interface 
from  the  layer  under  consideration  is  taken  as  zero  while  the  fluid  of 
the  layer  is  assumed  to  be  Incompressible  and  inviscid  and  the  flow 
irrotatlonal.  Moreover,  we  again  define  the  equilibrium  configuration 
of  a  liquid  layer  as  stable  whenever  any  initial  disturbance  of  suffic¬ 
iently  small  amplitudes  and  velocities  leads  to  a  motion  in  which  the 
an^jlltudes  and  velocities  remain  arbitrarily  small.  Also  we  assume  that 
the  eigenfunctions  corresponding  to  the  free  vibration  problem  are  com¬ 
plete,  which  allows  us  to  express  the  motion  of  the  fluid  due  to  any 
initial  disturbance  as  a  linear  combination  of  the  eigenvibrations.  We 
therefore  may  deduce  as  the  necessary  and  sufficient  stability  condition 
that  all  the  eigenvalues  (squares  of  the  eigenfrequencies)  be  real  and 
positive. 

For  layers  of  uniform  thickness  spread  over  various  portions  of  the 
exterior  or  interior  of  cylinders  and  spheres  the  stability  criteria 
could  be  presented  in  terms  of  inequalities  involving  the  magnitude  of 
either  the  unidirectional  body  force,  g,  or  the  nondlmensional  stability 
parameter 
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in  which  T  denotes  the  surface  tension,  p  the  fluid  density  and  a 
the  radius  of  the  equilibrivim  interface  surface.  With  these  layers  it 
is  consistently  the  first  antisymmetric  eigenvibration  which  leads  to 
the  stability  criterion,  since  the  corresponding  elgenfrequency  ceases 
first  to  be  real  when  we  vary  the  body  force  (or  stability  parameter) 
and  start  out  with  a  value  for  which  the  equilibrium  is  stable.  We 
anticipate  that  this  fact  also  prevails  in  layers  with  a  variable  thick¬ 
ness  and  therefore  limit  ourselves  to  determining  only  the  elgenfrequency 
of  the  first  antisymmetric  eigenvibration. 

While  exact  analytical  solutions  to  the  small  vibration  problem  of 

layers  with  a  nonuniform  thickness  can  in  general  not  be  obtained,  it  is 

nevertheless  possible  to  determine  rigorously  one  single  eigenvibration 

(12) 

in  certain  cases  by  means  of  a  simple  semi-inverse  method.  ^ 

choosing  to  find  the  first  antisymmetric  free  oscillation  with  the  aid 
of  this  method  and  comparing  the  results  for  various  layers,  we  can 
study  the  Influence  of  the  thickness  variation  on  the  stability. 
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II.  BASIC  EQUATIONS 


We  consider  layers  of  an  incompressible  liquid  which  are  bounded  by 
a  rigid  wetted  surface  and  an  interface  surface  S^^.  Ihe  interface 
separates  the  fluid  of  the  layer  under  consideration  from  a  neigh¬ 
boring  gas  or  fluid  of  zero  density.  Any  motion  of  the  liquid  is  assumed 
as  Irrotational,  which  allows  the  velocity  tT  to  be  expressed  with  the 
aid  of  a  velocity  potential  ♦  as 

r=  -V®  (l) 

Denoting  the  Laplacian  by  the  equation  of  continuity  has  the  form 


=  0  (2) 

We  consider  the  presence  of  a  body  force  k"  per  unit  mass  which  can  be 
derived  from  a  potential  Q  such  that 

r=  -vn  (3) 


Assuming  the  fluid  as  invlscid  and  having  a  density  p,  a  pressure  p, 
we  may  express  the  equilibrium  in  the  sense  of  d'Alembert  in  the  form 


'  ^  1 

p  +  p  1 

Integrating  Eq.  (4)  and  denoting  by  F(t)  an  arbitrary  function  of  time 
we  obtain  Euler's  equation 


p+p(n-||+|v^)  =  F(t)  (5) 

Without  loss  of  generality  we  may  Include  F(t)  in  the  time  variance  of 

In  addition  to  this  we  may  linearize  (5)  for  small  flow  velocities 

2 

by  neglecting  the  term  (l/2)pv  ,  which  reduces  equation  (5)  to 


The  kinematic  boundary  condition  of  the  fluid  layer  can  be  written  in 
the  form 
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m. 


0 


(7) 


where 

S  =  0  (8) 

represents  the  equation  of  the  bounding  surface  in  question. 

Denoting  the  surface  tension  by  T  and  the  external  pressure  by 
p  ,  the  dynamic  equilibrium  of  the  interface  surface  S.  can  be 
expressed  as 


s  2TH 


(9) 


where  and  R2  are  the  principal  radii  of  curvature  and  H  the 
mean  curvature.  (The  radii  of  curvatvire  are  taken  as  positive  if  the 
center  of  curvature  lies  within  the  fluid. ) 
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III.  SEMI-INVERSE  METHOD  OP  FINDING  EXACT  EIGENVIBRATIONS 


In  investigating  the  problem  of  the  free  oscillations  of  an 

incompressible  fluid  in  an  axially  symmetric  container,  considering 

exclusively  flat  interface  surfaces  and  neglecting  the  effects  of 

(12) 

surface  tension,  B.  A.  Troesch'  '  made  use  of  a  semi-inverse  method 
vhich  yields  a  rigorous  solution  for  one  of  the  eigenvibrations.  We 
apply  the  basic  idea  of  this  method  to  our  fluid  layers  in  the  following 
manner- 

Instead  of  prescribing  the  geometry  of  a  layer  completely,  we 
merely  select  the  shape  of  the  Interface  surface  in  a  suitable  manner. 
For  appropriately  chosen  solutions  to  Laplace's  equation  which  contain 
sufficiently  many  free  parameters  it  is  possible  to  satisfy  the  inter¬ 
face  boundary  conditions  in  the  presence  of  surface  tension  and  body 
forces.  By  establishing  the  streamlines  defined  by  such  solutions  and 
requiring  the  solid  supporting  surface  to  constitute  a  streamline,  we 
are  in  a  position  to  review  rigorous  vibration  solutions  corresponding 
to  a  class  of  liquid  layers  with  the  same  Interface  and  different  solid 
supporting  surfaces  and  thus  different  thickness  distributions. 

It  should  be  noted  that  the  selection  of  the  shape  of  the  interface 
together  with  the  body  forces  and  surface  tension  define  a  certain 
ambient  pressure  distribution,  the  feasibility  of  vhich  can  be  dealt 
with  separately- 
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IV.  EIGENVIHEIATIONS  AND  STABILITY  OF  TWO-DIMENSIONAL  LAYERS  WITH  FLAT 
INTERFACE 

1.  Rigorous  Solutions  to  Free  Vibration  Problem 

We  refer  the  Layer  to  Cartesian  coordinates  (x,  z)  such  that  for 
static  equilibrium  the  Interface  surface  Is  characterized  by  z  =  0 
(See  Fig.  l).  During  the  motion  resulting  from  a  small  disturbance  the 
Interface  surface  S^  deviates  from  Its  equilibrium  position  by  a  small 
distance  t)  and  the  equation  for  assumes  the  form 

S^(x,  z,  t)  =  z  -  ^(x,  t)  =  0  (10) 

With  Eq.  (lO)  the  mean  curvature  H  of  the  Interface  surface  Is  given 
by 


We  assume  that  ^(x,  t)  as  well  as  Its  derivatives  are  of  first  order 
small  and  neglect  all  terms  which  are  of  second  or  hl^er  order  small. 
Thus  Eq.  (ll)  reduces  to 

2H  =  -  (12) 

bx 

We  shall  restrict  ourselves  for  the  present  to  unidirectional  body  forces 
with  the  potential 

D  =  gz  (13) 

Substituting  Eqs.  (12),  (13)  and  (6)  Into  the  dynamic  boundary  condition 
(9)  we  obtain 

P  (x,  z,  t)  -  gzj  -  Pg(x)  =  -T  ^  (ll^) 


-  6  - 


With  Eqs.  (l)  and  (lO)  the  kinematic  boundary  condition  (?)  can  be 
written  in  linearized  form  as 

(W) 

By  differentiating  (l4)  with  respect  to  time  and  making  use  of  (15) 
we  can  combine  the  kinematic  and  dynamic  boundary  conditions  at  the 
interface  surface  into  a  relation  for  the  velocity  potential 


P 


= 


With  the  aid  of  Eq.  (2)  we  can  write  this  combined  boundary  condition 
also  in  the  form 


P 


5^4  S4^ 


(16) 


Neglecting  consistently  terms  which  are  of  second  and  higher  order  small, 
we  may  require  (l6)  to  be  satisfied  at  the  equilibrium  interface  z  =  0, 
rather  than  at  the  disturbed  interface  z  =  ^(x,  t). 

Since  we  are  seeking  solutions  in  the  form  of  harmonic  oscillations 
we  separate  the  space  and  time  variables  according  to 

4(x,  z,  t)  =  t(x,  z)e^"^‘^^  (17) 


and  thereby  reduce  the  initial-boundary  value  to  an  eigenvalue  problem 
with  ^(x,  z)  as  eigenfunction  and  a  as  eigenvalue.  By  substituting 
(17)  into  (2)  and  (I6)  we  find 


and 


-of 


+  T 


=  0 


Jz=0 


(IB) 


(19) 
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We  note  that  the  eigenvalue  problem  for  t|r  Is  of  a  special  type  since 
the  eigenvalue  appears  exclusively  in  the  boundary  condition  and  not 
in  the  differential  equation  and  since  the  boundary  condition  (19)  con¬ 
tains  derivatives  of  higher  order  than  does  the  differential  equation. 
With  the  transformations 


X 


1 


z  = 


z_ 

h 


where  h  is  a  characteristic  length  of  the  layer,  we  can  rewrite 
Eq[s.  (18)  and  (19)  in  terms  of  nondlmenslonal  space  parameters 


8x2  = 


(IS') 


oz  35 


=  0 


(19') 


z=0 


where 


P  =  and 

®  P6h2 


(20) 


For  the  sake  of  convenience  we  replace  again  x  by  x  and  z  by  z 
with  the  understanding  that  x  and  z  are  now  nondlmenslonal  length 
parameters.  The  eigenvalue  problem  can  thus  be  presented  in  the  form 


8x2  ^^2= 

=  0 

z=0 


“  rs 

oz 


(18) 

(19") 


We  consider  now  the  class  of  exact  solutions  for  t  which  can  be  obtained 
by  separating  the  space  variables: 

f(x,  z)  =  X(x)-Z(z)  (21) 

By  substituting  (2l)  into  (I8)  we  deduce 
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£ 

X 


(22) 


p 

whereby  k  may  assume  positive  as  well  as  negative  values.  We  consider 
first  k  0.  According  to  (22)  t  must  be  of  the  form 

t(x,z)  =  [A  sinh(kx)  +  B  coeh(kx)]  [C  sin(kz)  +  cos(k*)]  (23) 


If  k  0.  By  requiring  expression  (23)  to  satisfy  the  interface  bovuadary 
condition  (19")  we  obtain 


-3  +  k(l  -  MJt^)C  =  0  (2U) 

Since  the  solid  supporting  surface  of  the  layer  must  constitute  a 

streamline,  Is  given  hy 

Vt*n  =  0  (25) 


or 


df 


dz 


bf 


dx  =  0 


(26) 


where  tJ  Is  the  normal  vector  of  the  surface.  Integrating  (26)  by  making 
use  of  Eqs.  (2l),  (22)  and  (25),  we  find 


S2(x,z)  =  X'(x)-Z'(z)  -  E 

=  k^[A  cosh(kx)  +  B  8inh(kx)]-[C  cos  (kz)  -  sin(kz)  ]  -  E  =  0  (27) 


where  E  is  a  constant. 

If  we  restrict  oxirselves  to  layers  which  are  symmetric  with  respect 
to  the  z-axis,  we  have 

B  =  0 

Hence  the  eigenfunction  and  the  solid  surface  Sg  are  given  by 

t(x,z)  =  A  sinh(kx)  [C  sin(kz ) +  coB(kz) ]  (28) 

and 

S2(x,z)  =  co8h(kx)[C  C08(kz )  -  Bln(kz) ]  -  E  -  0  (29) 
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where  E  is  a  new  constant.  We  specifically  consider  the  following 
three  cases : 


(a)  E  =  0 

In  this  case  the  layer  is  of  uniform  depth  and  bounded  by  the  solid 
surfaces  x  =  ±Xq  and  z  =  -a  as  illustrated  in  Fig.  2.  Ihe 
corresponding  value  for  k  is  purely  imaginary 

k  =  li  (30) 

and 

Xq  =  ~2^'~  £  >  0,  n  =  integer  (3l) 

C  =  -1  tanh(ia)  (32) 

Substituting  (30)  and  (32)  into  (24)  and  (28)  we  obtain  a  simple 
relation  for  the  eigenvalue  a  and  the  corresponding  eigenfunction: 

P  =  —  =  /(l  +  |ii^)  tanh(ia)  (33) 

O 

t(x,  z)  =  D  sln(ix)  coshr/(a  +  z)]  (34) 


(b)  E  /  0 

The  parameters  C  and  E  in  Eqs.  (29)  may  be  determined  for  ex- 
anqple  by  requiring  the  surface  to  pass  through  the  points 
(0,  -a)  and  (xq,  O).  Thus  we  find  for  an  arbitrary  k  (real  or 
imaginary ) 

P=-^  =  k(i-:ik^)  ^  ■  (35) 

g  \  r-  /  coBh(kXQ)  -  co8(ka) 

Eco8h(kxQ)  -  cos(ka)  ”j 

os(kz) - -  8ln(kz)  -co8h(kxQ)  =  0 

♦(x,z)  =  D  8lnh(kx)  [cos  k(a  +  z)  -  coB(kz)co8h(kXQ) ]  (36) 

Representative  layers  corresponding  to  various  values  of  k  are 
shown  in  Fig.  3. 
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(c)  k  =  0 


In  this  case  the  solution  to  (l8)  is  of  the  form 


t(x,  z)  =  (Ax  +  B)(Cz  +  D)  (37) 

and  the  Interface  boundary  condition  ( 19” )  is  satisfied  if 

C  =  3D  (38) 


Considering  only  layers  which  are  symmetric  with  resi>ect  to  x  =  0, 
we  require 

B  =  0  (39) 

Ihe  integration  of  (26)  leads  to  the  expression 


S2(x,z) 


Choosing  the  values  of  (d/c)  and  E  such  that  the  surface  Sg 
contains  the  points  (0,  -a),  (xq,  O)  we  obtain 


C 

D 


ihl) 


and  layers  as  illustrated  in  Fig.  h. 


k,2.  Stability  Criteria 

Assuming  that  the  liquid  layer  under  consideration  has  finite 
dimensions  and  a  complete  set  of  eigenfunctions  * 

we  may  express  the  motion  of  the  fluid  due  to  any  Initial  dlst\irbance 
in  the  form 


from  which  we  immediately  derive 

o^^real,  o^  >  0,  n  =  1,  2,  ... 
as  the  necessary  and  sufficient  stability  criterion. 


(42) 


(43) 
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We  apply  this  criterion  to  the  eigenvibration  obtained  in  a  rigorous 
manner  in  the  previous  section  and  ve  consider  again  the  three  different 
cases:  E  =  0,  E  ^  0,  and  k  ^  0. 

(a)  E  =  0 

Requiring  o  >  0  in  (33),  we  deduce 


g(l  +  ni^)  >  0 

(44) 

With  Eqs.  (20),  (3l)  and 

b 

*0  "  h 

(45) 

Where  2b  is  the  width  of  the  layer  [See  Fig.  2]  we  can  reduce 
(U4)  to  the  classical  result 


Pg 

T 


> 


(b)  E  0 

From  Eq.  (3$)  follows 

o  >  0  if  gk(l-  ►ik^)sin(ka)  >  0  (47) 

We  note  that  Eqs.  (35)  and  (36)  do  not  change  if  we  replace  k  by 
-k  and  may  therefore  restrict  ourselves  to  positive  values  of  k. 
Since  the  slope  of  the  solid  surface  Sg  at  x  =  x^,  z  =  0  is 
given  by 


dz _ Blnh(4xo)Blrihfka) 

dx  ~  coBh(kxQ)[co8h(kxQ)  -co8(ka)] 

z=0 

its  sign  is  the  same  as  that  of  8in(ka). 

Ibis  means  that  we  have  a  convex  solid  surface  for  8ln(ka)<  0 
and  a  concave  one  for  8in(ka)>0.  Hence  considering  positive  as 
well  as  negative  values  for  g  we  may  restrict  ourselves  to 
Bin(ka)>0  without  loss  of  generality.  We  obtain  as  stability 
condition 
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a  >  0 


If 


^  ^  [sin(ka)  >  0] 


Similarly,  we  deduce  for  k  =  ii 


0  >  0 


(c)  k  =  0 

In  this  case  it  immediately  follows  from  Eq.  (li-l)  that 

a  >  0  if  g  >  0  (50) 

For  a  real  k  the  layers  are  only  stable  if  the  body  forces  exceed  a 
certain  positive  limit  value  while  for  an  imaginary  k  the  bodj-  forces 
may  come  close  to  a  certain  negative  limit  value  before  instability  is 
reached.  We  conclude  from  this  that  the  variation  of  the  thickness  of 
a  liquid  layer  may  significantly  influence  the  stability  criterion. 
Besides  this,  we  may  utilize  this  fact  in  designing  containers  of  fluids 
for  near -zero  gravity  conditions  in  such  a  manner  that  the  fluid  adheres 
to  a  certain  area  of  the  container. 
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V.  EIGENVIBRATIONS  AND  STABILITY  CRITiEIA  OF  IHREE-DIMENSIONAL  LAYERS 
WITH  FLAT  INTERFACE 


As  illustrated  in  Fig.  5,  we  refer  the  layer  tT  cylinder  coordinates 
z,  r,  6  such  that  for  equilibrium  the  Interface  surface  S^^  is  defined 
by  z  =  0.  Denoting  the  deviation  of  the  interface  from  its  equilibrium 
position  by  ^{r,  Q,  t)  we  may  write  for  the  surface 

S^(z,  r,  6,  t)  =  z  -  5(r,  0,  t)  =  0  (5l) 
and  its  mean  curvature  in  linearized  form 


As  before,  we  assume  the  presence  of  body  forces  whose  potential  is  de¬ 
fined  by  Eq.  (13).  The  combination  of  the  dynamic  and  kinematic  boundary 
conditions  leads  to 


We  separate  the  time  and  space  variables  according  to 
®(z,  r,  0,  t)  =  r,  0)e^'^  ^ 


(54) 


If  we  refer  the  length  parameters  again  to  a  characteristic  length  h 
we  may  interpret  z  and  r  as  nondimensional  quantities  and  define  the 
eigenvalue  problem  for  f  by 


^2^  1  1  ^2^  32^ 


=  0 


z*0 


Exact  solutions  for  f  can  be  obtained  in  the  form 


(55) 

(56) 
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♦(z,  r,  6)  =  F(z,  r)e 


Imd 


(57) 


which  requires  P  to  satisfy  the  equations 


.  1  ap  . 

p"*'? 


"•«-S 


=  0 


(58) 

(59) 


z=0 


Considering  first  such  functions  which  permit  an  additional  separation 
of  variables,  we  find 

F(z,  r)  =  A[cosh(kz)  +  B  sirih(kz)]J  (kr)  (60) 

m 

where  J  (kr)  denotes  the  m'*'^  order  Bessel  function  of  the  first 
m 

kind  and  k  the  separation  parameter.  The  substitution  of  (60)  into 
(59)  leads  to 

P  =  =  k(l  +  Mk^)B 


or 


a  =  k(l  +  nk^)B 


The  stability  conditions  are 


(61) 


k>0,  B^O, 

h 

and  the  corresponding  equation  for  Sg 

2  f  J  iiar) 

k  >  0,  m  >  0,  in[sinh(kz)  +  B  cosh(kz) ]  -  k  J  dr+  C  =0 

■|r 


(62) 


(63) 


Layers  with  solid  surfaces  defined  by  (63)  are  illustrated  in  Figure  6. 
We  consider  the  following  special  cases 
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(a)  =  0,  =  0 

Requiring  Sg  to  contain  the  points  (0,  -a)  and  (r^,  0)  we 
have  in  this  case  a  layer  of  uniform  thickness  bounded  by  the  svirface 
z  =  -a  and  the  cylinder  r  =  r^  (See  Fig.  ?)•  Ibe  relation  for  the 
eigenfrequencies  corresponding  to  mode  shapes  as  defined  by  (60)  is 

a  =  ^  k(l  +  iik^)tanh(ka)  (64) 


where  the  possible  values  of  k  are  given  by 


J^(hro)  =0  m  =  1,  2,  ... 


(65) 


From  (64)  we  find  the  stability  condition 


T  ^  ^ 


(66) 


which  indicates  that  the  lowest  value  for  k  is  most  significant :  m  =  1, 

krQ  =  1.85. 

1^}  ■  2 

In  this  case  the  velocity  potential  can  be  written  in  the  form 


f(rj  Z;  9)  =  A  r'’^  (z  +  B)e^™®,  m  0 


(67) 


The  interface  boundary  conditon  (59)  yields 


(68) 


and  the  equation  of  the  solid  surface  Sg  becomes 


/ 1  2  \  12 
Sg  =  m  z^  +  Bzj  -  i  r  +  C  =  0 


(69) 


By  requiring  the  surface  Sg  to  pass  through  the  points  (r  =  0, 
z  =  -a)  and  (r  =  r^,  z  =  O)  we  find 


B  = 


2^2 
a  m  +  r^ 

2ma 


(70) 
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and 


2^2 

2  ®  ®  ^0  2  2 

S„  =  mz  + - ^  z  +  -  r  =  0  (71) 

«=  a  ^ 

From  (68)  and  (70)  we  obtain  ae  stability  condition 

g  >  0  (72) 

which  holds  Independently  of  m,  Tq  and  a.  Representative  layers 
defined  by  (71)  are  shown  in  Pig.  8. 

As  in  the  two-dimensional  case,  we  conclude  that  the  variation  of 
the  thickness  can  markedly  Influence  the  stability  of  three-dimensional 
layers  with  a  flat  Interface,  particularly  for  body  forces  in  the  near¬ 
zero  region.  We  also  note  that  for  some  solid  supporting  surfaces  Sg 
(container  walls)  layers  of  this  type  are  stable  if  the  body  forces  are 
larger  than  a  certain  negative  limit  value  while  for  other  surfaces  Sg 
the  body  forces  have  to  be  positive  for  stability,  ^toreover,  we  may 
again  suggest  the  feasibility  of  exploiting  such  a  behavior  in  the 
design  of  fluid  containers  for  near-zero  gravity  conditions. 
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VI.  ON  THE  STABILITY  OF  LAYERS  WITH  UNIFORMLY  CURVED  INTERFACE 


6.1.  TVo-Dlmenslonal  Layers  with  Curved  Interface 

We  refer  the  layer  to  cylindrical  coordinates  such  that  the  equil¬ 
ibrium  configuration  is  characterized  by  r  =  a  as  shown  in  Figure  9* 

In  the  presence  of  a  disturbance  the  equation  for  the  interface  surface 
may  be  written  in  the  form 

S^(r,  0,  t)  =  r  -  [a  +  ^(0,  t)]  =  0  (73) 

In  terms  of  cylindrical  coordinates,  the  uniaxial  body  forces  are  now 
defined  by  the  potential 

n  =  gr  cos  0  (7^) 

and  the  linearized  mean  curvature  given  by 

■  I  ^  i  ^  <■'5) 


By  combining  the  dynamic  and  kinematic  boundary  conditions  as  before,  we 
obtain 


6$ 

.1.  er  -  U 

T 

at"  *  57  «■= 

2 

a 

r=a 

_  . 

We  again  separate  the  time  and  space  variables  according  to 
4>(r,0,t)  =  *(r,0)e^'^^ 


(76) 


(77) 


and  may  interpret  r  as  nondimenslonal  length  parameter  by  using  "a" 
as  the  reference  length.  Thus  we  arrive  at  the  eigenvalue  problem  for 
defined  by 


—  .  1  3^  .  1_  ^ 

^j.2  r  5r  ^2  ^g2 


0 


(78) 
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(79) 


=  0 


r=l 


where  3  and  n  now  abbreviate 


\x  =  and  P  = 

Pga^  8 


(80) 


The  streamlines  can  be  obtained  by  integrating  the  relation 


^  ^  =  7  Se  ^ 


(81) 


Any  function  of  the  form 
N 

Ir.-' 

n^l 


Xi 

,e)  =  ^  (r  ”  +  r  ’^)[Bjj  sinCij^e)  +  cos(in©)] 


(82) 


with  arbitrary  in  aad  N 

we  may  choose  i  =  n  and 
n 

3 

Hr,e)  =  ^  (r"  +  A^ 
n=l 

We  require  =  0,  n  =  1,2,3  and  thus  restrict  our  attention  again  to 
the  antisymmetric  modes  of  oscillation.  The  substitution  of 

3 

t(r,0)  =  y*  A^r''^]sln(n0) 

iSi 

into  the  Interface  boundary  condition  (79)  leads  to  the  relations 


represents  a  solution  to  (78).  In  particular, 
N  =  3: 

r  '^)tB^  8in(n0)  +  cos(n0)]  (83) 
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V®2  '  ’■ 

Ag  =  1  -  23\ 

V®2  =  ^  +  2p[l  -  + 

Vl/®2  =  -  2P[1  -  X.(P 

^3  =  1 

in  which  \  Is  an  arbitrary  parameter.  If  we  select 

V  -  1 

"  P  + 


6n)]  \ 

+  6|i)] 


(84) 


(85) 


we  limit  our  analysis  to  a  family  of  layers  for  which  we  can  readily 
derive  a  stability  condition  by  replacing  the  abbreviations  P  and  n 
in  (85)  by  their  original  form  and  requiring  a  to  be  positive.  We 
deduce 


a  >  0, 


a  <  0, 


(86) 


With  (85)  the  differential  equation  for  the  stream  lines  becomes 

dr  r(\(r -r'^)sln  0-t-2[r^-  (l- 2\p)r'^]sln  26  +  3\(r^- r~^)sln  59) 
\(r  +  r  ^)cos  0  +  2[r^+  (1 -2\p)r  ^]coB  20+ 3\(r^+ r^)cos  30 


Representative  layers  with  such  streamlines  are  illustrated  in  Figure  10. 

6.2.  Two-Dimensional  Layers  with  Slightly  Curved  Interface 
For  such  layers  we  may  write 

cos  0  »»  1  (88) 

and  as  interface  condition 
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Ihls  condition  Is  satisfied  by 


f(r,e)  =  B(r^  +  Ar‘^)sin  iG 
for  arbitrary  £  if 

A  =  i[l  +  ^i(i^  -  1)]  -  g 
i[l  +  -  l)]  +  p 


(90) 


(91) 


The  corresponding  stream  lines  are  given  by 
B[r'^  -  Ar"'^ ]cos{i0)  -  C  =  0 


(92) 


where  C  denotes  a  constant.  %  choosing  C  =  0  we  arrive  at  a  family 
of  layers  with  uniform  thickness  and  the  same  stability  criteria  as  dis¬ 
cussed  in  Ref.  2.  For  C  ^  0  we  may  require  the  solid  boundary  to  contain 
two  prescribed  points,  e.g. 

9=0,  r  =  TqJ  0  =  ±a,  r  =  1 
which  determines  A  and  C/b  in  (92): 


A  = 


r^  -  cos(ia) 

-  cos(£a) 


(93) 


-£  £ 

r  -  T' 

C  0  0 


®  rZ^  -  co8(£a) 


cos(io:) 


(9^) 


By  comparing  (93)  and  (9l)  we  find  the  relation 


2£ 


a  =  i 

a 


1  +  -  2rQ  co8(^a) 


g  +  -2^  (/  -  1) 

_  pa  J 


(95) 


from  which  we  can  deduce  the  stability  conditions 


a  >  0  Tq  <  1:  >  1  -  / 

a  <  0  Tq  >  1:  2^  <1-1^ 

As  a  special  case  we  consider  i  =  1  and  A  =  0.  The  substitution  of 
(93)  and  (94)  into  (92)  yields  for  the  solid  boundary 

r  cos  0  =  Tq  (96) 

l.e.,  a  flat  sxirface.  The  corresponding  velocity  potential  is  given  by 
t  =  Br  sin  0  (97) 


Specific  examples  of  two-dimensional  layers  with  a  slightly  curved  inter¬ 
face  are  illustrated  in  Pig.  11. 


6.3.  Three-Dimensional  Layers  with  Curved  Interface 

In  terms  of  spherical  coordinates  we  may  define  the  disturbed 
Interface  by 

S^(r,  0,  cp,  t)  =  r  -  [a  +  £(0,  cp,  t)]  =  0  (98) 

where  r  =  a  represents  the  equilibrium  configuration  as  shown  in  Fig.  12. 
The  potential  of  uniaxial  body  forces  is  again  given  by  (7^).  For  small 
S  the  mean  curvature  may  be  written  in  the  form^^^^ 


2H  =  i 

a 


2  -  -  2f  + 
a 


^  cot  0  +  ^  +  — - 


_  ^  ^ 

sin20 


(99) 


A  combination  of  the  kinematic  and  dynamic  boundary  conditions  for 
the  Interface  leads  to 


^g  g  cos  0  )  -  %  (2 


+  cote  I  -  I  ^  „  I 

^  SrSe^  Sir!?0  arij^J 


=  0 


r*a 


(100) 
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The  physical  nature  of  the  problem  suggests  a  separation  of  the  variables 
according  to 


(101) 


In  which  m  denotes  an  Integer.  We  may  again  Interpret  r  as  the 
nondlmenslonal  length  parameter,  using  the  radius  "a"  of  the  equili¬ 
brium  shape  as  reference  length.  The  eigenvalue  problem  for  ♦  Is  now 
defined  by 


(102) 


(103) 


The  nondlmenslonal  quantities  p  and  n  are  again  defined  by  (80)  and 
the  stream  lines  by  (8l).  Exact  solutions  to  the  eigenvalue  problems 
are  for  example 


A  [r^*^ 
n 


+  (cos  e) 

i'n 


(104) 


In  which  represents  the  associated  Legendre  function  of  degree  m 

and  order  ’  The  order  p^  Is  arbitrary  while  the  degree  m  Is 

Integer . 

We  consider  a  particular  family  of  layers  to  which  (104)  represents 
an  exact  solution  by  selecting 

m  =  1  ,  p^  =  1,  2 
or 

*  =  Aj^[(r  +  B^r"^)sln  0  +  +  B2r"^)sln  0  cos  0]  (105) 


This  fixnctlon  satisfies  the  Interface  condition  If 
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(106) 


-1,  B2  =  §, 

From  (106)  and  (80)  follows 


i. 

5P 


=  -2-  fi 

5A2  a 


For  Ag  >  0,  a  >  0  we  find  as  the  stability  condition 

g  >  0 

layers  representative  for  this  family  are  shown  in  Pig.  13. 


(107) 


6.4.  Three-Dimensional  Layers  with  Slightly  Curved  Interface 
For  such  layers  we  have 

cos  0  **  1,  sin  6^0  (108) 

and  therefore  as  the  continuity  equation  and  the  interface  condition 

2 


(109) 


-  II  .  ^ 


2 

+r  — 7!  +  r  — T 


br 


dr" 


(no) 


r=l 


Solutions  to  Eq.  (109)  which  satisfy  the  Interface  condition  (llO)  can 
be  written  in  the  form 


t(r,0)  =  [rP  +  Ar-(P+l)]j  [n/p(p  +  l)'  9] 

m 

where 

-P +p{l+ U[p(p +1)- 2]} 

A  =  i.i—  -  -  1.  «  I  — 

+P  +(p  +  l)(l  +  n[p(p+l)-  2]) 


(111) 


(112) 


and 


p  >  0  . 


The  streamlines  corresponding  to  (ill)  are 
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In 


^p(p+l70 

.(P*!)  .  *  (EU)  .-Pj  .f  ^ 

0 


dX  +  C  (m  >  O) 


(113) 


with  C  denoting  a  constant. 

For  C  =  0  we  have  a  family  of  three-dimensional  layers  of  uniform 
thickness  as  Illustrated  in  Fig.  l^t.  ?y  defining  the  solid  surface  as 


r  =  r^j  0  =  *  a  (a  small) 


(114) 


and  requiring 

37  =  50  =  ° 

we  obtain 

p  +  1  0 

and 


J‘  [n/p(p  +  l)  a]  =  0 
m 

The  comparison  of  (ll5)  with  (ll2)  leads  to 
1  P(l-r^^^) 


(115) 


(116) 


a  =  - 
a 


1  + 

p  +  1  0 


g  +  —  [p(p  +  l)  -  2] 
pa 


(117) 


with  rQ  <  1,  and  p  >  0,  we  have 


g  +  — 2  +  l)  -  2]  >  0 

pa 

and  therefore  as  stability  condition 


g  >  [2  -  p(p  +  l)] 


(118) 


pa 


Similarly,  we  find  for  >  1 
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(119) 


g  <  [2-  p(p  +l)]  -~ 

It  should  be  noted  that  for  a  given  Ot  the  parameter  p  is  determined 
by  (ll6).  According  to  (ll8)  and  (II9)  we  obtain  the  stability  criterion 
with  the  smallest  value  of  p-  The  latter  is  given  by  the  smallest  root 
of 

[s/p(p  +  i)'  a]  =  0  (120) 

The  physical  interpretation  of  this  fact  is  that  among  all  eigenvibrations 
the  first  antisymmetric  oscillation  leads  to  the  stability  criterion  since 
the  corresponding  frequency  approaches  zero  first  in  varying  the  parameters 
T,  p  and  g. 

Hie  stability  criterion  for  an  axially  symmetric  layer  with  a  flat 
Interface,  a  uniform  thickness  and  a  radius  b  as  shown  in  Fig.  7  is 


g  >  -13.7 


T 

p(2b)^ 


while  that  of  the  sUgJitly  cxirved  layer  of  uniform  thickness  and  the  same 
interface  area  (see  Fig.  l4)  is 


g  > 


15.7~|  .  T 

(2bfJ  P 


As  in  the  two-dimensional  case'  we  observe  that  the  curvature  has  a 
destabilizing  effect. 


If  C  ^  0  in  Eq.  (113),  we  can  interpret  (ll2)  as  the  frequency 
equation  for  arbitrary  values  of  p  and  A: 


{' 


g  +  [P(P  +1)  -  2] 


pa 


'} 


For  stability  we  deduce 


(121) 


(p  t.  >  0 

1  +  A 


a  ^  0  , 


g  ^  [2  -  p(p  +  1)1 

pa 
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a  ^  0  , 


(122) 


<c. 


&  ^  — 5  [2  -  p(p  +  l)] 

pa 


Representative  examples  of  slightly  curved  three-dimensional  layers  are 
Illustrated  In  Fig.  15* 

If  we  assume  p(p  1)8  «  1>  we  may  obtain  approximate  solutions  to 
the  streamline  Eq. (113)  In  the  form 


J  3^  ^  -  S&TT7_ 


2  -I 


For  m  =  1  we  have 


/ 


®  J  (X) 
/  m 


tn 


1  -  2l£ 


and  thus 


S,  =  [r^^-  -  A  r-p]  [l  .  Sllfii  9=]^  - 


C  =  0 


(123) 


Requiring  the  solid  surface  to  contain  the  points  (r  =  1,  0  =  0()  and 


(r  =  Tq,  0=0)  we  find 


A  “ 


(124) 


and 


Sg  “  I  r 


.  1^1. ^  .  (,-P.  ^1)  I"!.  0 


(125) 


For  p  =  1  and  A  =  0  we  have  an  approximation  to  the  Important  case 
of  a  layer  on  a  flat  surface 
r  cos  0  =  tq 
as  shown  In  Figure  15b. 
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VII.  CONCLUSION 


The  results  of  this  study  clearly  lnd.icate  a  strong  dependence  of 
the  stability  of  all  layers  (two-  and  three-dimensional)  on  the  thickness 
variation,  especially  in  the  realm  of  near-zero  gravity  conditions.  Of 
particular  significance  is  the  fact  that  for  certain  solid  supporting 
surfaces  (container  walls)  the  fluid  layers  are  stable  only  if  the  uni¬ 
directional  body  force  is  either  positive  or  exceeds  a  certain  positive 
value,  while  for  other  surfaces  stability  prevails  for  limited  negative 
values  of  the  body  force.  taking  advantage  of  this  fact  in  the  design 
of  containers  of  fluid  for  near-zero  gravity  conditions,  it  appears  pos¬ 
sible  to  enforce  the  permanent  adherence  of  the  fluid  to  a  certain  area 
of  the  container. 

It  should  be  noted  that  we  have  not  Introduced  any  contact  angle 
considerations  In  explicit  form  in  our  analysis,  and  therefore  we  have 
to  require  that  the  dimensions  of  the  interface  surface  be  large  compared 
with  the  effective  range  of  the  contact  forces.  Also,  it  should  be  point¬ 
ed  out  that  our  stability  criteria  are  conservative  only  if  the  "appropri¬ 
ately  chosen"  solution  to  Laplace's  equation  actually  leads  to  the 
fundamental  eigenvi brat ion*  whose  eigenfrequency  ceases  first  to  be  real 
in  varying  the  body  force  or  stability  parameter.  Partial  reassurance  in 
this  regard  is  offered  by  the  fact  that  our  results  include  as  special 
cases  the  known  stability  criteria  for  layers  of  uniform  thickness  with 
either  flat  or  uniformly  curved  interface  surfaces. 


Vibration  with  the  lowest  eigenfrequency. 
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Fig.  I.  Reference  Frame  for  Two-Dimensional 
Layer  with  Flat  Interface. 


z 


Fig.  2  Two-Dimensional  Layer  whose  Free 

Oscillations  have  the  Vslocity  Potential 
♦(x,z,t) «  Dsin(/x)cosh[/(a+z)]*e‘'^* 
where  /•[(2n-l)/2]r(h/b). 
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-K|  S„(a-I,  Xo-0.7,  k»2i) 

(b)  k  imaginary 


Fig.  3  Two-Dimensional  Layers  with  Free  Oscillations 
Defined  by  ^(x,z,t)*Dsinh(kx)[cosk(a<«-z)-‘COs(kz) 
co8h(kxo)]e'^*;  (a)  k  real,  (b)  k  imaginary. 
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Fig.  5  Reference  Frame  for  Three-Dimensional 
Layer  with  Fiat  Interface. 
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Fig.  6  Three-Dimensional  Layers  with  Free 
Oscillations  Defined  by 

<I>(r,z,0,t)  =  A[cosh(kz)  +  Bsinh(kz)]j|(kr)e 


Fig.  7  Three-Dimensional  Lover  whose  Free 
Oscillations  can  be  characterized  by 

^(r,z,0,t)  -  Acosh(ka+kz)  Jni(kr)e'^"'®'*‘‘^^' 

where  Jm(kro)sO 
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Fig.  8  Three-Dimensional  Layers  with  Free 
Oscillations  Defined  bv 
•Wr.  z,9,t) = Ar[z+(a'fro*)/2a]e*  * 


Fig,  9.  Reference  Frame  for  Two-Dimensional 
Layer  with  Curved  Interface. 
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e^o 


Fig.  10  Two-Dimensional  Layers  with  Oscillations 
Defined  by 

4>(r,0,t)  =  B{x(r+|)sin6+[r*+(l-2iSX)r‘*]sin2^+ 
X(r*+r"*)sin3^}e'^*  where  X* 
{l/I(<7-o/g  +6(T^go»)]} 
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e^o 


(b) /  =  l  ro=cosa 
<t=Brsin0e'*^^ 

Fig.  I  la,b.  Two-Dimensional  Layers  with  Slightly 
Curved  Interface. 
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0=0 


(d)  <I>  =  B[(^/+(^)^]cos/0e'‘^’  where  /  =  y  or 

<I>  =  B[(y +(]^f^]sin/0e*^*  where 

Fig.  lie, d.  Two-Dimensional  Layers  with  Slightly 
Curved  Interface 


Fig.  12.  Reference  Frame  for  Three-Dimensional 
Layers  with  Curved  Interface. 


Fig.  13.  Three-Dimensional  Layers  with  Oscillations 
Defined  by 

^  s  A[(r-r"*)sin^+Ag(r*+%r"*)sin0co88]e'^^*^^^' 
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Fig.  14  Three-Dimensional  Layers  with  Slightly 
Curved  Interface  and  Uniform  Thickness. 
The  Fundamental  Oscillation  of  such 
Layers  is  Characterized  by 

where  a)  •  0 
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e«o 


A“ 


and  >/ptp+l)  8  « I 


^S2»(p*1»ro*cosa) 
(b)  <t  =  BrJ,(y28)e'‘*-"'^*' 


0=0 


Fig.  15.  Three-Dimensional  Layers  with  Slightly 
Curved  Interface  and  Oscillations  Defined 
as  above. 
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